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1: Introduction 

The ubiquitous Brownian motion remains an outstanding paradigm in modern physics. 
Some representative, but by no means an exhaustive hst of general references ("founding 
papers", reviews and apphcations) are presented in [1-19]. Here we present the Langevin 
formulation for a Brownian carrier in uniform and static external fields. Some recent work 
on charged Brownian particles is referenced in [20-47]. In our previous work on this mat- 
ter, our approach hinged on the resolution of Kramers and/or Smoluchowski equations 
[24,26,28,45,46,47], and recently we began to tackle Langevin 's formulation of this prob- 
lem [47]. Here we explore the latter, in order to study the relaxation of the Brownian 
carrier towards a steady state, given electrical and magnetic external static and uniform 
fields. In section 2 we present the solution of Langevin 's equation including the above 
mentioned fields. In section 3 we present our results for the nonequilbrium temperature 
relaxation to the "hot" steady state temperature (as modified by the electric and magnetic 
fields). The computed final "hot" regime temperature is compared to long time existing 
results ([48] with no magnetic field present) and with our previous results, with the mag- 
netic field contribution, via Kramers and Smoluchowski equations [26,28,45]. In section 4 we 
present an application, namely a simple yet relevant Brownian model (with no adjustable 
parameters) for GaAs carrier mobility [49-56]. The multivalley band structure, and the "hot" 
carrier steady state temperature obtained in the previous section are the essential ingredi- 
ents for the appearance of a negative differential conductivity regime, in good quantitative 
agreement with well known experimental results. Furthermore our model incorporates the 
magnetic field contribution hitherto not considered. Finally , in section 6 we present our 
concluding remarks and outline further work. 

2: Langevin equation for a Brownian charged particle 

We briefiy present the Langevin formalism for a free charged Brownian particle [10-12,18,19], 
with mass m, and charge q immersed in an homogeneous thermal reservoir at temperature 
Tr. It is essentially Newton's equation for the particle with two contributing forces: the 



first, a systematic dissipative force Stokes like (linear in the particle's velocity) and the 
second a rapidly fluctuating random force, 

rj\r rn 

m— = Fs + F^ = -7V + F'-(t) r = - (1) 
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The formal solution is 

v(t) = exp Li"! v" + ^ / dti exp (^^) F'-(ti) (2) 

with initial condition v° = v(0) and r the collision time. The random force has solely 
statistical properties: zero average and white noise correlations , given by the averages 

(F^(t)> = {Fl{h)F]{h)) = 2-kBTnS,,6{h - h) (3) 

where the correlation strength is such that the asymptotic average kinetic energy satisfies 
the equipartition theorem, in thermal equilibrium with the thermal reservoir(fluctuation 
dissipation theorem), and given by 

hrn {w\t ^ oo)> = heTR = ^mVS (4) 

Following [28,45,46,47] (and with a slightly different notation) we now consider the Brown- 
ian carrier (charged particle) under the influence of homogeneous external, time independent, 
electric and magnetic fields; the electric contribution is given by Feiec = Q'E and the magnetic 
contribution (Lorentz's velocity dependent force) Fmag=^Q'v x B. The total external force 
is given by 

F(v)= Feicc+Fmag(v) =qE-mco XV CO = — B (5) 

Let us define a tensorial Stokes force by adding the Lorentz contribution to the usual 
Stokes force, as 



TS = 



-7V— -mw X V = —mA "^v (6) 



where u is the usual cyclotron frequency, the magneto mobility tensor is M =m~^A with A 
a tensorial collision time, that can be cast into the form (when operating over an arbitrary 
vector V) 



A(t,uW = T 



V+rV X uj+t'^u (w ■ V) 



l + T^OO 



2,. ,2 
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In particular notice the familiar form for the case B =B^ 



A(t,u;) 



1 + T^U 
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By defining such a tensorial Stokes force, Langevin's equation now reads 



m- 



Itt 



-mA"V+gE + F'(t) 



with formal solution [57-60] 



(8) 



(9) 



v(t) = exp(-A-H)v° + A(l-exp(-A-^t))— + — f rftiexp (A-^(ti - t)) F'(ti) (10) 

f^ "^ Jo 

Using Cayley-Hamilton theorem, and Putzer [59] and Apostol [60] results, after a lengthily 
but straight forward calculation we obtain 



exp 



exp (A 4) = ao(t) + ai(t)A ^ + a2(t)A ^ 
t\ f . 1 ,. , ,.o,. . 1 
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(11) 
(12) 



3: Evolution of the Effective (nonequilibrium) temperature towards a hot 

steady state. 



Let us define the nonequilibrium effective temperature Tef(t) as 

hsT^,{t) = lm{^r\t)) (13) 

and consider "thermal like" initial velocity conditions, i.e. configurational averages 0^, over 
initial velocities are given by Maxwell like distributions 

(v°>, = lm{v^v';)^=hnnS,, (14) 

where Tq < Tr (Tq > Tr) defines cold (hot) initial distributions. Furthermore we define 
dimensionless electric and magnetic fields e and b as 






h =TUJ eh=ebcos9 u = \uj\ (16) 

and 6 the angle between fields. Furthermore we introduce a dimensionless time in colli- 
sion time units and denoted hereafter as t = t/r. Finally, by considering equation (3), 
some trivial integration and basic vector and matrix algebra, the expression for the effective 
temperature is given by 



Tef(t) = To exp (-2t) +Tn{l- exp (-2t)) + T^e' ('1±^1^ + T^{t) + T^it) ] (17) 
with 



2exp(-t) ff_^^ , 26 ^,^^,\^,o^ , _2 



ri(t) = "^ ' ' I I cos6t + - T^sinbt ] sin^^ + cos^^ ) (18) 
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r2(t) = exp(-2t) ( 1 + 6^cos^6'H 2 cos 6t sin 6t sin^ 6* + cos^ 6* ) (19) 

(1 + 6^) 
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The Brownian carrier gas evolves towards a stationary state with a nonequihbrium tem- 
perature greater that the reservoir temperature, thereby the name hot carrier. 



«-^"i)r"-^n^Ti;^J^ (20) 

This result includes the magnetic field effect. In the Kramers Smoluchowski scheme [26,28,45] 
we derived the expression 

Tf^"'{t ^ oo) 1 (l+ffm£l\ ( V«(x) V 

where a long tail (1/t non exponential tail, see also [25])and the spatial inhomogeneity of 
the carrier's distribution corrects the electric field contribution. Both the long tail and the 
density diffusive effect smears out at longer times rendering both expressions equivalent. We 
believe equations (20) and (21) represent novel results at least in the Brownian context, as 
far as the magnetic field contribution is concerned. Shockley obtained an expression for the 
hot carrier temperature [48-50] for null magnetic field and in agreement with our result in 
equation (20) for 6 = 0. The absence of the electric field renders the equilibrium (reservoir) 
temperature, regardless of the magnetic field value. For non zero electric field values, the 
magnetic field modulates the hot carrier's temperature with maximum value for parallel 
fields (independent of the magnetic field value) and minimum value for perpendicular fields. 
Figures 1 to 5 we plot the effective temperature evolution, with unit time given by the 
collision time constant r. Solid lines represent the zero magnetic field and dashed lines the 
magnetic and angle value on display in each figure. The reservoir temperature is fixed at 
400°K and two values are chosen for Tq namely 200°K (cold) and 600°K (hot) initial velocity 
distributions. The cold (hot) cases are the lower (upper) curves starting at t = 0. We observe 
that after a few collision time units the carrier's temperature reaches the stationary value 
that decreases as the magnetic field deviates from the parallel configuration, similarly for 
the non parallel case as the magnetic field increases. The transient effects are damped 
oscillations, a larger effect for intermediate magnetic field values. The dimensionless field 



values where chosen from typical material parameters for GaAs (see next section), rendering 
e = 1 corresponds to 1000 volts/cm and 6=1 corresponds to one Tesla. 
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Figure 1: Effective temperature versus time, electric field e = 3.0. Solid lines represent the 
zero magnetic field case 6 = 0, dashed lines the case 6 = 0.1, ^ = 0. Hot initial condition 
given by Te/(0) = 600^ K and cold initial condition by Te/(0) = 200^ K . See main text for 
time and field units. 
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Figure 2: Effective temperature versus time, electric field e = 3.0. Solid lines represent the 
zero magnetic field case b = 0, dashed lines the case b = 10.0, ^ = 0. Hot initial condition 
given by Te/(0) = 600^ K and cold initial condition by Tef{0) = 200^ K . See main text for 
time and field units. 
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Figure 3: Effective temperature versus time, electric field e = 3.0. Solid lines represent the 
zero magnetic field case 6 = 0, dashed lines the case 6 = 0.1, 9 = 7r/2. Hot initial condition 
given by Tef{0) = 600^ K and cold initial condition by Tef{0) = 200^ K . See main text for 
time and field units. 
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Figure 4: Effective temperature versus time, electric field e = 3.0. Solid lines represent the 
zero magnetic field case 6 = 0, dashed lines the case b = 2.0, 6 = 7i/2. Hot initial condition 
given by Te/(0) = 600°i^ and cold initial condition by Tef{0) = 200°i^. 5*66 main text for 
time and field units. 
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Figure 5: Effective temperature versus time, electric field e = 3.0. Solid lines represent the 
zero magnetic field case 6 = 0, dashed lines the case h = 10.0, 6 = 7i/2. Hot initial condition 
given by Te/(0) = 600°i^ and cold initial condition by Tef{0) = 200°i^. 5*66 main text for 
time and field units. 
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4: Simple Brownian carrier model for GaAs. Negative differential 

conductivity. 

Gallium arsenide (GaAs) is a compound of the elements Gallium and Arsenic. It is a 
III/V semiconductor, and is used in the manufacture of devices such as microwave frequency 
integrated circuits, monolithic microwave integrated circuits, infrared light-emitting diodes, 
laser diodes, solar cells and optical windows. The band structure consists of a muhivalley 
landscape. Based on data from references [49-56] we model the electron mobility /i (propor- 
tional to the drift velocity Vo), combining the high mobility F valley with low mobihty L 
valley (six satellite valleys per F type valley). The carrier mobility is obtained via standard 
canonical average restricted to these two type of valleys, 

, = tE^^ (22) 

-rr + -TL 

with the usual expressions as given in the literature [49-56] 



jWa = /^oa/— ^t; Pa^ ga'ml^'^expi-—^] a = F,L (23) 



where /io is proportional to a typical mean velocity equation (10) and with relevant typical 
parameters given by: gi^ = 6gr, ttli^ ~ lOmp and A = Ei^ — E^ — 0.3eV. At room temperature 
300°i^ we have A/kBTf> ^ 12. The drift velocity (mobility times electric field) is computed 
in the electric field direction as (in arbitrary units). 

V, = iil(ii2£)5 (24) 

With typical values for GaAs we compute the last equation using equations (10) and (22) 
yielding the material parameter free equations 



_ 1 / l + 60exp(-12e-^) \ 
^^-7ilvl + 13Oexp(-120-i)J^'"' ^25) 

1 + b^cos^g 
^ = i + b^ (26) 



= 1 + OMQe^ (27) 
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notice that for zero gap we obtain a typical Caughey Thomas [56] expression for the mobihty. 

,.„ ^ o.47n 

^ Vl + 0.04fie2 ^ ^ 

Caughey Thomas mobihty modehng is a non hnear phenomenological fitting procedure. 
Our general result is linear in the electric field (in the linear regime, the current is a 
transport coefficient times the field), incorporates the magnetic field, and the nonlinear 
electric field dependence is due to the non equilibrium temperature intrinsic field dependence. 
As mentioned in the previous section, the dimensionless field values where chosen from 
typical material parameters for GaAs, rendering: e = 1 ^ E = 1000 volts/cm and b = 1 ^ 
B = ITesla. 

In Figures 6 to 9 we plot the drift velocity in arbitrary units versus electric field. The 
solid upper line is the zero magnetic field case and the lower solid lines is the correspond- 
ing Caughey Thomas case (A = 0). The dashed lines are corresponds to the angle and 
magnetic field values on display in each figure, the dashed lines maximum decreases as the 
magnetic field value increases. For all cases presented we find a region of negative differential 
conductivity (Gunn effect) 

dVDJe) ^Q 
de 

for electric field larger than the critical value Cc ~ 4000 Volts/cm, ( — ^^^ ~ '-' ) ^-s it is 
well known for this compound [49,50,52-54]. From Figures 6-9 and considering equations 
(25-27), as the magnetic field value is increased we notice that edb) increases and Voiedb)) 
decreases while the effective temperature 0(ec(6), 6, 6) decreases. This pattern becomes more 
pronounced as we move from parallel fields towards perpendicular fields. The available ex- 
perimental data for a related compound [55], where the negative differential conductivity 
region was probed with a magnetic field, seems to corroborate our findings, in a very qual- 
itative fashion. In the higher electric field regions, say e > 10 where our results deviates 
from the experimental data, the effective temperature O is very large, quite probably other 
scattering mechanisms should be incorporated, rendering a more involved temperature de- 
pendence for the mobility (see equation (23)), and more than one collision time constant 
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should be considered. We include very large magnetic field values solely to probe the pattern 
described above. 
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Figure 6: Drijt velocity versus electric field. Upper solid line for zero magnetic filed (or 
arbitrary magnetic filed with 9 = 0). Dashed lines in decreasing order for b = 0.5, 1.0 and 
2.0, with 9 = tt/4. Lower solid line case of null gap and zero magnetic field. See main text 
for velocity and field units. 
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Figure 7: Drift velocity versus electric field. Upper solid line for zero magnetic filed (or 
arbitrary magnetic filed with 9 = 0). Dashed lines in decreasing order for b = 5.0, 20.0 and 
50.0 (the last two cases overlap on the scale used) , with 9 = 7r/4. Lower solid line case of 
null gap and zero magnetic field. See main text for velocity and field units. 
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Figure 8: Upper solid line for zero magnetic filed (or arbitrary magnetic filed with 9 = 0). 
Dashed lines in decreasing order for b = 0.5, 1.0 and 2.0 , with 9 = 7i/2. Lower solid line 
case of null gap and zero magnetic field. See main text for velocity and field units. 
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Figure 9: Drift velocity versus electric field. Upper solid line for zero magnetic filed (or 
arbitrary magnetic filed with 9 = 0). Dashed lines in decreasing order for b = 5.0, 20.0 and 
50.0 (the last case almost indistinguishable from the e axis, on the scale used) , with 9 = 7r/2. 
Lower solid line case of null gap and zero magnetic field. See main text for velocity and field 
units. 
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5: Concluding remarks 

We presented the Langevin formulation for a Brownian carrier under uniform and static 
external electric and magnetic fields. From the solution to the associated Langevin equation, 
we computed the relaxation of the carrier's effective (nonequilibrium) temperature towards a 
(hot) steady state regime with a nonequilbrium field dependent temperature. The latter was 
compared with well known existing results. We believe our present result in the Langevin as 
well as our previous results in the Kramers Smoluchowski scheme, incorporate the effect of 
the magnetic field hitherto not considered. Then we presented a simple yet relevant Brownian 
model to account for the negative differential conductivity behavior on the GaAs compound, 
again incorporating the magnetic field effects hitherto not considered. Discussions of results 
and Figures are presented at the end of sections 3 and 4. 

Our future work includes the incorporation of diffusive effects on the effective temperature, 
as discussed in section 3, the inclusion of other scattering mechanisms into the mobility as 
discussed in section 4, and incorporate within the Langevin Formalism chemical reactions 
[61-62] and photovoltaic effects [63] as discussed for example in [45] within the Kramers 
Smoluchowski context. 

Acknowledgments: We thank Profs. Ricardo Paupitz and Mariano de Souza for valuable 
comments. This work was partially supported by CNPq (Brasil). 
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